Low-energy fixed points of random Heisenberg models 
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The effect of quenched disorder on the low-energy and low-temperature properties of various two- 
and three-dimensional Heisenberg models is studied by a numerical strong disorder renormalization 
group method. For strong enough disorder we have identified two relevant fixed points, in which the 
gap exponent, w, describing the low-energy tail of the gap distribution, P(A) ~ A" is independent 
of disorder, the strength of couplings and the value of the spin. The dynamical behavior of non- 
frustrated random antiferromagnetic models is controlled by a singlet-like fixed point, whereas for 
frustrated models the fixed point corresponds to a large spin formation and the gap exponent is 
given by u ~ 0. Another type of universality classes is observed at quantum critical points and 
in dimerized phases but no infinite randomness behavior is found, in contrast to one-dimensional 
models. 

PACS numbers: 05.50.+q, 64.60.Ak, 68.35.Rh 
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I. INTRODUCTION 

The Heisenberg model plays a central role in the 
theory of magnetic ordering^ and the two-dimensional 
(2d) antiferromagnetic (AF) model has been intensively 
studied motivated by its relation to high-temperature 
superconductivity^. According to the Mermin- Wagner 
theorem-" 7 , no long-range order (LRO) can persist at fi- 
nite temperatures in the homogeneous Heisenberg model 
if d < 2. At zero temperature, the LRO of the clas- 
sical ground state is reduced by quantum fluctuations. 
This effect is particularly strong in (quasi)- Id AF models 
and gives rise to the complete destruction of Neel-type 
LRO. Fluctuations enhanced by quenched randomness 
and frustration can further destabilize LRO, resulting 
in disordered ground states even in higher-dimensional 
systems. In various experiments, in which quasi-two- 
dimensional magnetic materials that can appropriately 
be described by the 2d HAF were diluted with static non- 
magnetic impurities (Mg or Zn in La2Cu04, and Mg in 
K2C0F4 or K2MnF4), a disorder-induced transition from 
Neel order to a spin liquid was observed: If the impurity 
concentration is larger than a critical value the LRO is 
destroyed 3 ^. 

The behavior of Heisenberg antiferromagnets (HAF-s) 
in the presence of quenched randomness is generally very 
complex and present understanding of this is not com- 
plete. Most of the theoretical results have been obtained 
for Id models, many of them by a strong disorder renor- 
malization group (SDRG) method introduced originally 
by Ma, Dasgupta and Hu for the random S — 1/2 AF 
spin chain£. Fisher— has shown that the SDRG method 
leads to asymptotically exact results in the vicinity of a 
quantum critical point, which corresponds to the chain 
without dimerization. At the quantum critical point, the 
ground state can be described by the notion of a ran- 
dom singlet (RS) phase, which consists of effective sin- 



glets of pairs of spins that are arbitrarily far from each 
other. Fisher's SDRG treatment has been extended to 
the dimerized phases that turned out to be equivalent 
to quantum Griffiths phases^. The SDRG method has 
also been applied for random S = 1— and S = 3 /2m spin 
chains and for various random spin ladder models 1 ^. In 
general, the Haldane gapped phases stay gapped for weak 
disorder, while they become gapless and often form RS 
phases for strong disorder. 

To study the singular properties of the S — 1/2 Heisen- 
berg model with mixed ferromagnetic (F) and AF cou- 
plings, the SDRG method has to be modified. In Id, the 
presence of ferromagnetic couplings leads to the forma- 
tion of large spin clusters in the RG treatment, with an 
effective moment that grows without limits as the energy 
scale is lowerediii As a consequence, the ground state 
properties of random Heisenberg chains with mixed AF 
and F coupling and of those with only AF couplings are 
different. The presence of large effective spins in the low- 
energy limit was also observed for random AF spin lad- 
ders with site dilution. 12 

Not many theoretical investigation of the effect of 
quenched disorder in higher-dimensional random HAF-s 
exist, and those that have been done are almost exclu- 
sively restricted to dilution on the square lattice. Quan- 
tum Monte Carlo studies of the HAF on a diluted square 
lattice show that LRO disappears at the classical per- 
colation poinljii While in earlier investigations a novel, 
5-dependent critical behavior was foundi 3 -, recent stud- 
ies identify the transition as an S-independent classical 
percolation transition with the well known exponents 14 . 
Another work studied the ± J Heisenberg (quantum) spin 
glass and found that for a concentration of F bonds p > 
p c sa 0.11 the Neel-type LRO in the ground state vanishes 
and is replaced by a so-called spin glass phase 15 . Within 
the spin-glass phase, the average ground state spin, S^ot, 
scales as Stot ~ VN, and the gap as AE ~ 1/N, where 
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N is the number of spins>i& 

In this paper we study the effect of randomness 
in higher dimensional HAF-s by means of the SDRG 
method. In particular, we consider the low-energy behav- 
ior of frustrated and non- frustrated systems in 2d and 3d. 
As we will mention in the next section the pure (i.e. non- 
random) versions of these models have a ground state 
that has either AF or dimer LRO or is disordered, i.e. in 
a spin-liquid state. By calculating the gap distribution 
and cluster formation within the SDRG scheme we shall 
characterize the change of the pure ground state struc- 
ture of the pure systems by the effect of the disorder. 

The paper is organized as follows: The models and 
their phase diagrams for non-random couplings are pre- 
sented in Sec. II. The SDRG method and its different 
low-energy fixed-points for (quasi)-l(i systems are dis- 
cussed in Sec. III. Results of the SDRG method on dif- 
ferent 2d and 3d models are presented in Sec. IV and 
discussed in Sec. V. 



II. THE MODELS AND THEIR PHASE 
DIAGRAM FOR NONRANDOM COUPLINGS 

We start with the Hamiltonian of a nearest-neighbor 
spin-1/2 AF Heisenberg model: 

Ht = JS ^' . (!) 

{kk') nn 

where J > and the summation runs over nearest neigh- 
bor(nn) pairs, (k,k'), of a regular lattice. In Id, Neel- 
type LRO is destroyed by quantum fluctuations and the 
system with half-integer spin value S show quasi-long- 
range-order (QLRO), i.e. correlations in the ground state 
decay algebraically 1819 . A similar behavior can be ob- 
served in AF spin ladders with an odd number of legs 20 . 
Both systems have a gapless excitation spectrum and in 
finite chains of length L the gap vanishes algebraically 
with a dynamical critical exponent z c = 1, which is char- 
acteristic for a quantum critical point (QCP): 

AE CI (L) ~ 1/L'' . (2) 

Quantum fluctuations play a different role in AF spin 
chains with an integer spiniSt and for spin ladders with 
an even number of legs 2 ^. These systems show a topolog- 
ical string order, which is accompanied by exponentially 
decaying correlations and by a finite gap in the spectrum. 

One can approach a 2d geometry by successively in- 
creasing the number of legs of AF spin ladders. The 
resulting square lattice has a qualitatively different low- 
energy behavior: The effect of quantum fluctuations is 
weaker and the ground state shows Neel-type LRO 22 . 
Compared with the classical ground state the sub-lattice 
magnetization for S = 1/2 is reduced by about 40%. 
Generally in ordered AF phases the excitation spectrum 
is gapless. In finite e?-dimensional system of linear size 
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FIG. 1: The J\ — J2 model (a) and the dimerized model (b) 
on the square lattice. 

L - according to spin-wave theory and analysis of the 
non-linear sigma model - the gap behaves as^: 

AE od (L) ~ 1/L d . (3) 

Frustration generally leads to a further reduction of the 
Neel LRO. Frustration of geometrical origin is present 
in the triangular lattice, where the sub-lattice magne- 
tization is about 50% of its classical valued. In more 
loosely packed frustrated lattices, such as in the kagome 
latticed and in the square lattice with crosses2&(see, how- 
ever Ref. I22}) or in the 3d pyrochlore lattice 28 , the LRO 
completely disappears and the systems have a disordered 
ground state. The correlations are short ranged and one 
finds a finite triplet gap in which a continuum of singlet 
excitations exists. In the case of the kagome lattice these 
extend down to the ground stated. 

Competing interactions are another source of frustra- 
tion which can also lead to disordered ground states. As 
an example, we consider the AF J\ — Ji model with first- 
( Ji) and second-neighbor ( J2) interactions, described by 
the Hamiltonian: 

H = H 1 + H 2 , (4) 

where 

{kk') nnn 

and the coupling in H\ (Eq.QJ) is denoted as J = Ji (see 
Fig. [TJi) . In 2d, there are at least three phases as shown 
in Fig 2. a. For small frustration, J2/J1 = p, the system 
possesses AF LRO, whereas for large frustration the sys- 
tem goes to the collinear state, in which ferromagnetically 
ordered columns of spins are arranged antiferromagnet- 
ically. In the range 0.34 < p < 0.60, the ground state 
is disordered and the spectrum is gapped for all types 
of excitations^. According to recent numerical studies^ 
there are probably several quantum phases in this region, 
separated by different type of quantum phase transitions. 

Finally, we introduce a dimerization into model Q. 
We consider the square lattice, denote a lattice site k by 
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FIG. 2: Phase diagrams of square lattice HAF models, a) For 
the J1-J2 model with varying frustration, p = J2/J1, there are 
three regions: the ordered AF phase and the ordered collinear 
(CL) phase, separated by a disordered spin-liquid (SL) region, 
b) In the dimerized model, the AF and dimer (D) ordered 
phases are separated by a quantum critical point at a c . 

its two coordinates, k = (i, j) and define 

-ffdim = — ^ JOiS2i,jS2i+l } j ■ (6) 

i,j 

The dimerized model is then described by the Hamilto- 
nian H = Hi + Hdim and has a layered structure, see 
in Fig. ^d. Its phase diagram is shown in Fig. 03 as a 
function of the dimerization parameter < a < 1. For 
a < a c = 0.686 the ground state has AF LRO, whereas 
for a > a c the system is in an ordered dimerized phase, 
in which spin-spin correlations along vertical lines ap- 
proach different limits if the distance between the spins 
is odd or even lattice sites, respectively. In the dimer- 
ized phase, there is a finite gap which vanishes at a c as 
AE ~ (a~a c y with v = 0.71, characteristic for the uni- 
versality class of the 3d classical Heisenberg model2L22. 
We note that dimerization with another topology has 
been studied recently in Ref. 13^. 

The random Heisenberg models we investigate in this 
paper include the 2- /3-dimensional AF model on the reg- 
ular lattice Q), the dimerized AF model © in 2c?, geo- 
metrically frustrated AF models on the triangular lattice 
as well as on the kagome lattice, the 2-/3- dimensional 
Ji — J2 model and the 3d AF-F models. We are interested 
in how the phase diagrams in Fig- [3 will be modified due 
to the presence of strong quenched randomness. 

III. THE SDRG METHOD AND ITS 
LOW-ENERGY FIXED POINTS IN Id MODELS 

The basic ingredient of the SDRG method in Heisen- 
berg models is a successive decrease of the energy scale of 
excitations via a successive decimation of couplings. We 
start with a S = 1/2 HAF model in which the strongest 
coupling is, say J23, the one between lattice sites 2 and 
3 (c.f. Fig. I3J) - If J23 is much larger than its neighbor- 
ing couplings J12, </i3, J24 and J34, the spins at 2 and 3 
form an effective singlet and are decimated. The effective 



coupling between the remaining sites, 1 and 4 in second 
order perturbation theory is given by: 

iff = A (Jl2 - Jl3 P 4 ~ J24) , A(S = 1/2) = 1/2. 
J23 

(7) 

In a chain geometry the couplings J13 and J24 would not 
be present and the resulting RG flow always generates 
AF couplings. However, for extended, not strictly Id 
objects, some of the generated new couplings can be fer- 
romagnetic (e.g. if J12 < J13 and J34 > J24 or vice versa) 
and therefore the decimation rules have to be extended. 
If at one RG step an F bond turns out to be the strongest 
one, its decimation will lead to an effective spin S = 1. 
In the following steps, the system will renormalize to a 
set of effective spins of different magnitude interacting 
via F and/or AF couplings. 

For higher dimensional systems, the basic decimation 
processes are the singlet formation in Eq.Q) and the ef- 
fective spin (cluster) formation. To specify the latter, let 
us consider three spins Si , S2 and S3 with interactions 
fulfilling IJ23I ^S> I «7i2 1 , j «/x3 1 • In the action of the RG, 
the two original spins £2 and S3 form a new effective 
spin of magnitude S — |$2 ± #3] representing the to- 
tal spin of the ground state in the two-spin Hamiltonian 
H23 = J23S2S2, where the positive (negative) sign refers 
to an F (AF) coupling. The corresponding energy gap, 
A, between the ground state and the first excited state 
in the Hamiltonian H23 is given by A = | J23 1 + S3) 
and A = J23(|5 2 — 5 3 | + 1), for an F and AF coupling, 
respectively. If J23 > (AF) and S2 = S3, it follows an 
effective singlet formation as described above. If 5 ^ 0, 
within first order perturbation theory the new coupling 
between Si and S23 is given by 

J eS = C12 Jl2 + Cl 3 Jl 3 , (8) 

with 

S(S + 1) + S 2 (S 2 + 1) - S 3 (S 3 + 1) 

C12 = =rr= 

25(5+1) 

and 

5(5 + 1) + 5 3 (5 3 + 1) - 5 2 (5 2 + 1) 

C13 = =— = . 

25(5+1) 
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FIG. 3: Singlet formation and decimation for a spin config- 
uration that does not have a chain topology and typically 
occurs in higher dimensional systems. 
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At each RG step, we find the pair of the spins with 
the largest energy gap A that sets the energy scale, fi, 
and decimate them according to renormalization rules 
described in J7J) or JSJ. A detailed derivation of these 
renormalization rules can be found in Ref. l34l 

The fixed point of the RG transformation for lattices 
that do not have a chain geometry may depend on their 
topology, the original distribution of bonds, the strength 
of the disorder, etc. We briefly summarize the existing 
results for spin chains and ladders since it might be help- 
ful for analyzing the RG results in higher dimensional 
systems. 

In the case of the random AF chain (which does nei- 
ther have F bonds nor frustration), the RG procedure 
described above runs into an infinite randomness fixed 
point (IRFP) corresponding to a random singlet phase. 
In this phase the renormalized clusters are singlets, thus 
the total magnetic moment is zero, and the energy and 
length scales are related via 



In ft - L 



1/2 



(9) 



which means that the dynamical exponent is formally 
infinite. 

A dimerized S = 1/2 chain with random AF even (J ) 
and odd (J ) couplings shows dimer order and the low- 
energy behavior is controlled by a random dimer fixed 
point at which the dynamical exponent, z, is finite and 
a continuously varying function of the strength of the 
dimerization measured by 5dim = [In J c ] a v — [In J ]a-^^- 
At this fixed point, the low-energy-tail of the distribution 
of the effective couplings, J e , is given by: 



p{j c ,n)dj c 



1 (Je 



1 1 <u 



~TT , (10) 



for 5(j; m > 0. This random dimer phase is a Griffiths 
phased and we refer to it as a Griffiths fixed point (GFP) . 
At this GFP, the gap of finite chains of length L obey a 
distribution similar to Ea. H10|) : 



P L (A) = L Z P(L Z A) ~ £2(1+") A a 



(11) 



which is characterized by the gap exponent, to. As a 
consequence of Eg. (|llfl - which holds in any dimension, 
several dynamical quantities at a GFP are singular and 
the characteristic exponents can all expressed via uj. For 
example the susceptibility %, the specific heat C v (at a 
small temperatures T), and the magnetization m (in a 
small field h), behave as: 



m(h) 



•OJ+1 



(12) 



In the Griffiths phase there is a simple relation between 
the dynamical exponent, z, and the gap exponent, uj, 
which can be obtained by the following phenomenologi- 
cal consideration 37 : If the Griffiths singularities are due 
to rare events (produced by the couplings) that give rise 
to localized low-energy excitations, the gap distribution 



should be proportional to the volume, 
From Ea. ((TT)> then follows: 



L" 



(13) 



which is consistent with the exact result in the random 
dimer phase in Ea. (|l(Jf> . However, if the low-energy exci- 
tations are extended the relation 113|) might not hold. 

In a spin chain with mixed F and AF couplings^, large 
effective spins, S e g, are formed at the fixed point of the 
transformation. The size of of these spin clusters scales 
with the fraction of surviving sites during decimation, 
l/N, as: 



S, 



off 



(14) 



The following random walk argument 11 gives C = 1/2: 
The total moment of a typical cluster of size N can be 
expressed as S e g = | SiLi where neighboring spins 

with F (AF) couplings enter the sum with the same (dif- 
ferent) sign. If the position of the F and AF bonds are 
uncorrelated and if their distribution is symmetrical, one 
has S cS oc N 1 / 2 , i.e. Eq.Q with ( = 1/2. 

A non-trivial relation constitutes the connection be- 
tween the energy scale f2 and the size of the effective 
spin: 



S'off ~ CI 



(15) 



where a numerical estimate of the exponent is k = 
0.22(1)11. Comparing Eq.lO with Eq.^, the relation 
between the length scale L ~ N 1 ^ (d = 1) and the en- 
ergy scale is: 



CI 



K 



1 

2k 



(16) 



where z is the dynamical exponent. The distribution of 
low-energy gaps, Pj,(A), has the same power-law form as 
in Eg. ( |11| >. Therefore from the scaling behavior of Pl(A) 
the gap exponent, u>, and the dynamical exponent, z, 
can be obtained. Due to the large moment formation 
the singularities of the dynamical quantities are different 
from those in the random dimer phase in Ea. (|12l) . i.e. at 
a GFP. Generalizing the reasoning in Ref. fill we obtain 
in d-dimensions: 

X(T) - T-\ C V (T) - T 2C( " +1) |lnT|, m{h) - fe^wS) , 

thus the singularities involve both exponents Q and uj. In 
the following, we refer to this type of fixed point as large 
spin fixed point (LSFP). 

AF spin ladders, although being quasi-one- 
dimensional, have a non-trivial, non-chain-like topology 
and during renormalization also F bonds can be gener- 
ated according to Eq.JJJ. Different random AF two-leg 
ladders were studied in Ref. with tho following 
results: If the disorder is strong enough the gapped 
phases of the non-random systems become gapless. The 
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low energy behavior is generally controlled by a GFP, 
where the dynamical exponent is finite and depends 
on the strength of the disorder. However, at random 
quantum critical points, separating phases with different 
topological or dimer order, the low-energy behavior is 
controlled by an IRFP. In diluted AF spin ladders also 
LSFP-s have been identified^. 

To close this section we summarize that in one- 
dimensional and in quasi-one-dimensional random 
Heisenberg systems there are two different types of low- 
energy fixed points, which are expected to be present in 
higher dimensional systems, too. Both for a GFP and 
for a LSFP, the low-energy excitations follow the same 
power-law form in Ea. (|ll|l from which the exponents, w 
and z can be deduced. At a GFP these two exponents 
are expected to be related through z = l(T^|) . On the 
other hand, for a LSFP, where the excitations are not 
localized, this relation probably does not hold. At such a 
LSFP there is a third independent exponent, £ involved 
in the dynamical singularities partially listed in Eg. (II 711 . 

In the next section we are going to study different two- 
and three-dimensional random Heisenberg models. In 
particular we will be interested in the possible difference 
in the low-energy fixed point for non-frustrated and frus- 
trated systems. Since extended (quasi-one-dimensional 
or higher-dimensional) random HAF models and Heisen- 
berg models with mixed F and AF bonds follow the 
same renormalization route, they could, in principle, be 
attracted by the same fixed points, but also new fixed 
points can emerge, as we will see. 



IV. RENORMALIZATION OF 
HIGHER-DIMENSIONAL SYSTEMS 

This section is the central part of our work, where we 
present our results for the ground state structure of vari- 
ous two- and three-dimensional random Heisenberg mod- 
els obtained by the numerical application of the SDRG. 
In practice we start with a finite system of linear size 
L with periodic boundary conditions and perform the 
decimation procedure up to the last effective spin (or 
decimate out the last spin-singlet). The energy scale cor- 
responding to the last decimation step is denoted by A. 
This procedure is performed for several thousand real- 
izations of the disorder and yields a histogram for A, 
which represents our estimate of the probability distri- 
bution Pi (A). From this we extract the gap exponent 
u> and the dynamical exponent z via the asymptotic re- 
lation given in Ea. l|ll|) . Moreover, from the average size 
of the effective spin at the last step, = [S c s]av, the 
cluster exponent, £, in Ea. (|14fl is deduced. The value of 
uj, z and £ will then be used to discriminate the different 
possible low-energy fixed points described in the previous 
section. 

Throughout this paper we use a power-law distribution 



for the random couplings < J < 1 for AF models: 

Pd(J) = Ij- 1 ^ , (18) 

where D 2 = [(lnJ) 2 ] av — [In J] 2 V denotes the strength 
of the disorder. Note that both the initial distribution 
of the couplings in Eq. l(T%|) and the final distribution of 
gaps in Ea. Hll|) follow power laws. If + 1) < D, the 
strength of disorder is reduced during renormalization, 
thus the low-energy random fixed point is a conventional 
one. More generally, for a conventional random fixed 
point ui > — 1. In contrast to this, at a IRFP the disor- 
der growths without limits, thus here formally u) = — 1 
and the dynamical exponent is infinite. We often use 
the uniform distribution, which corresponds to D = 1 in 
Ea. (|18() . For models with random F and AF couplings 
we take either a Gaussian distribution: 

P G (J) = -=L= exp(- J 2 /2a 2 ) , (19) 

or a rectangular distribution: 

P r {J) = 6(J- r + l/2)8(r + 1/2 - J) , (20) 

where 9(x) = 1, for x > and zero, otherwise. The latter 
distribution is symmetric for r = 0, whereas for r = 1/2 
we recover the uniform distribution of AF couplings in 
Ea.iP)) with D = 1. 

A. Two-dimensional models 

In the calculations for 2d we usually considered systems 
of linear size up to L — 32, but for some cases in which the 
convergence was faster we went only up to L = 10 — 16. 
The typical number of realizations were several hundred- 
thousands for the smaller sizes and several ten-thousands 
for larger systems for each value of D. At the first part 
we investigate non- frustrated models, such as the HAF 
on the square lattice with and without dimerization. In 
the second part of our study we consider frustration, the 
origin of which could be i) geometrical as for instance 
for the triangular and kagome lattice ii) due to a random 
mixture of F and AF couplings as for instance for the ± J 
spin glass model and iii) due to competition between first- 
and second-neighbor couplings as for the J\-Ji model. 

1. HAF on the square lattice 

We start with the renormalization of the HAF on the 
square lattice. The probability distribution of the gap 
calculated for a uniform bond distribution fEq. <|18[) with 
D = 1) is shown in Fig. 0] for different linear sizes. In a 
log-log plot the small gap region of the curve is linear, the 
slope of which, according to Eq. I|ll|) , corresponds tow + 
1. With increasing size one observes a slight broadening 
of the distributions indicating a decreasing effective gap 



6 



exponent which, however, seems to converge to a finite 
asymptotic value: 



u AF = 0.7(1), d = 2 



(21) 



During renormalization we observed simultaneously an 
effective singlet formation, thus in Eg. (|14l) one has £ = 0. 
Our estimate for the dynamical exponent satisfies the 
relation in Ea. (|13l) . yielding zaf = 1-2. Thus we con- 
clude that the low-energy fixed point of the system is a 
conventional, finite disorder Griffiths fixed point and the 
thermodynamical singularities are given by Eq. (|12fl . For 
other disorder strengths D we reach the same conclusions 
and our estimates for the gap exponents for each D agree 
with the value in Eci. (|21|) within the error bars. Thus 
the low-energy singular behavior of the 2d random HAF 
does not depend on the strength of disorder, in contrast 
to random quantum spin ladders^. 



2. Square-lattice HAF with dimerization 

Next we study the low-energy behavior of the dimer- 
ized HAF, as sketched in Fig. \%)p. For site and bond 
dilution the stability of the gapped, dimerized phase was 
recently investigated^. Here we consider the effect of 
strong AF bond disorder. In our calculation we used 
uniform initial randomness and performed the renormal- 
ization for several values of the dimerization parameter, 
a. The possible values of the two types of couplings were 
in the regions (0, 1] and (0, (1 — a)], respectively. For any 
value of a in the range < a < 1, we observed an effec- 
tive singlet formation and the estimated gap exponents u) 
and dynamical exponents z are found to satisfy the rela- 
tion in Eq. (|13fl . The extrapolated dynamical exponents 
as plotted in Fig.[5]seem to be approximately constant in 
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FIG. 5: Extrapolated dynamical exponent of the random 
dimerized HAF on the square lattice. The random AF and 
the random dimerized phases are separated by a cross-over 
region in which the dynamical exponent is minimal. 



two regions, which corresponds to the two phases of the 
pure model in Fig. For weaker dimerization the dy- 
namical exponent corresponds to the one of the random 
HAF, and for stronger dimerization z is approximately 
equal to the one of the disconnected two-leg ladder sys- 
tem, to which the case a = — 1 reduces, with z 1.07 10 . 
We expect that the dimer order is finite in the RD re- 
gion, whereas it is zero (or very small) in the random 
HAF region. Between the two regions, corresponding to 
the neighborhood of the phase transition point in the 
pure system in Fig. [2), the dynamical exponent drops 
to a minimal value. This cross-over could happen in a 
smooth, non-singular way, or in a sharp phase transition 
separating the random AF and the random dimer phases. 
Due to strong finite-size effects we could not discriminate 
the two scenarios. 

We note that z in the cross-over region behaves in 
the opposite way as in the dimerized ladders, where the 
dynamical exponent at the transition point in a finite 
system is maximal, and increases without limits^ for 
increasing system size, signaling an IRFP. In the two- 
dimensional case considered here the combined effect of 
critical fluctuations and quenched randomness seem to 
reduce the value of the dynamical exponent. Our cal- 
culations indicate that in the random dimer phase the 
low-energy behavior is controlled by a GFP and the dy- 
namical singularities are given by Ea. (|12fl . 



3. Randomly frustrated models (2d ) 



FIG. 4: Distribution of the energy gap of the square-lattice 
HAF with uniformly distributed random couplings, for linear 
sizes L = 8, 16, 24 and 32. The slope of the low-energy tail of 
the distributions is given by — (uj + 1) = —d/z. The straight 
line for L = 32 has a slope m —1.7. 



In this subsection we consider the Heisenberg model on 
the square lattice with a random mixture of F and AF 
couplings. This is a model for a quantum spin glasa^^ 
and we denote the corresponding fixed point as spin glass 
fixed point (SGFP), although we do not explicitely check 
for the existence of proper spin glass order in the ground 
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FIG. 6: Probability distribution of the energy gap on the 
square lattice with mixed F and AF bonds following a Gaus- 
sian distribution with a = 1. (The slope of the straight line 
is —1). Inset: Distribution of the spin moments. 



state (for instance via the calculation of the Edwards- 
Anderson spin glass order parameter—). As we will see, 
this fixed point differs from the other fixed points we 
found for non- frustrated models, so we feel that the use 
of this notation is justified. In particular we find a large 
spin formation proportional to L during RG procedure 
implying a ground state spin S oc V^V, which is remi- 
niscent of the spin glass behavior found ini£ii& for this 
model with alternative methods. 

First we report the results for the Gaussian random- 
ness in Ea. l|19|l . For this case the distributions of the gaps 
and of the effective spin moments are shown in Fig. 
The gap-distributions for different finite sizes have a very 
similar structure: they are merely shifted to each other 
by a constant proportional to IriL. The slope of the low- 
energy tail of the distributions is practically independent 
of the strength of disorder and in all cases the gap expo- 
nent is equal to: 



uj sg = 0, d = 2 , 



(22) 



within an accuracy of a few percent. From the finite- 
size scaling properties of the gap distribution, we infer 
that the relation in Ea. <|13|) is satisfied and therfore the 
excitations are localized, implying 



z SG = 2, d = 2 



(23) 



within an accuracy of a few percent. 

On the other hand, the distribution of the effective 
spin moments in the inset to Fig. shows a tendency to 
broaden with increasing system size and its average value 
has a linear L dependence, [/i/Jav ~ .42L. Therefore the 
moment exponent in Ea. ll4|) is 



Csg = 1/2, 



(24) 



We have repeated the above analysis using the sym- 
metric rectangular distribution in Ea. (|20|l both for the 



S = 1/2 and the 5 = 1 models and we obtained the 
same critical exponents as in the Gaussian case. Thus we 
can conclude that the low-energy behavior in randomly 
frustrated 2d models is controlled by the same SGFP, in- 
dependent of the type of randomness and the size of the 
spin. 



4- Geometrically frustrated models 

In this section we consider the HAF on two geometri- 
cally frustrated lattices that have qualitatively different 
ground states in the non-random case. The triangular 
lattice has finite AF long-range order and low-energy ex- 
citations behave as in Eq.©. In contrast to this, the 
ground state of the kagome lattice is disordered and the 
low-energy singlet excitations have a more complicated 
size dependence. 

We start with the HAF on the triangular lattice us- 
ing the power-law distribution in Ea. i|18|) for the ran- 
dom couplings. The distribution function of the gap is 
presented in Fig. [7| for different disorder strengths. The 
slope of the low-energy tail of the distributions is again, 
as for the randomly frustrated model of the last susec- 
tion, practically independent of the strength of disorder 
and in all cases the gap exponent is equal to u> — within 
an accuracy of a few percent. 

When calculating the moment of the spin clusters, we 
notice large spin formation during the action of the RG. 
From the size dependence of the average moment we ob- 
tain the exponent in Ea. H14|) £ = 1/2, independent of the 
strength of disorder. From the finite-size scaling proper- 
ties of the gap distribution, we infer that the relation in 
Eq. (|13|) is satisfied and therfore the excitations are local- 
ized, implying z$g — d = 2. Thus we can conclude that 



Triangular lattice, L=10 



o 




FIG. 7: Probability distribution of the energy gap for the 
triangular lattice HAF for different strength of randomness in 
Eq.flHJ- The low energy tail of the distributions, which has 
practically no finite-size dependence for L > 10, is consistent 
with the same gap exponent, u) = 0, implying a dynamical 
exponent zsa = d = 2. 



8 



Kagome, D=2.5, L=8 
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FIG. 8: Dynamical exponent of the random HAF on the 
dimerized kagome lattice with a randomness parameter D = 
2.5 calculated in finite systems having L 2 triangles, thus 3L 2 
sites. 



the thermodynamical quantities in the random triangular 
HAF obey the relations in Ea. lfT?)) . 

Next we focus on the kagome lattice and enlarge the 
parameter space by c onsidering the dimerized model, as 
introduced in Ref. y!J Couplings in up-pointing triangles 
(J) are different from those in down-pointing triangles 
(J'). Analyzing the results of the RG calculation as al- 
ready described for the triangular lattice, we obtain a set 
of gap, dynamical and moment exponents for different 
dimerization, 0.1 < J' / J < 1.5, and disorder strengths, 
D = 1, 2.5, and 5. In Fig. [S|we show our estimates for 
the dynamical exponents for D — 2.5, which are con- 
sistent with the SGFP result in Ea.l(22|). Also for other 
disorder strengths we find the same behavior and we con- 
clude that the low energy physics of the random kagome 
HAF is controlled by the SGFP and the thermodynamic 
singularities are described by Ea.Q17|). 



5. The J1-J2 model 

In our final example for the 2c? case, the source of frus- 
tration is the competition between first - J\- and second- 
neighbor - J 2 - couplings, which obey a power-law distri- 
bution in Eq. ifTHl) within the range of < J\ < J™ ax and 

< J2 < J™ 31 *, respectively. We have performed the pre- 
vious analysis at different points of the phase diagram, 
jmaxyjmax^ an( j f or different strength of disorder, D. In 
all cases we found that the relation in Ea. (|13fl is valid. 
As an illustration we show in Fig.[5]our estimates for the 
dynamical exponents for a disorder strength D = 5/3, 
which are consistent with the SGFP value in Ea. (|22|l in 
a wide range of 0.2 < J™ ax /J™ ax < 2.0. The same con- 
clusion holds for other disorder strengths in the range of 

1 < D < 5. During renormalization there is large spin 
formation and the calculated cluster exponent is consis- 
tent with = 1/2. Thus we can conclude that in the 



J f J 2 model, D=5/3,L=10 
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1.8 - 11 
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j 2 /j 1 

FIG. 9: Dynamical exponent of the J1-J2 model on the square 
lattice with a power-law randomness with D = 5/3. 

J\-Ji model the different phases in the pure model (AF 
and CL ordered, disordered SL) are washed out by strong 
disorder, and the whole frustrated region, J2/J1 > 0, is 
controlled by the SGFP. 

B. Three-dimensional models 

For the calculations in 3c? that we present now we con- 
sidered only systems of linear size L = 6,8, 10 and 12, in 
some cases we went up to L — 16. Larger system sizes 
were computationally not feasible. The typical number 
of realizations were several ten-thousand for each point. 
Due to the smaller system sizes the finite-size effects in 3c? 
are stronger than in 2c?. These finite size effects turned 
out to be too strong in the random HAF on the cubic 
lattice for a safe estimate for the gap exponent. We can, 
however, conclude that there is no large spin formation 
and the low-energy behavior is controlled by a conven- 
tional GFP. 



1. Randomly frustrated models (3d) 

We have studied models with mixed F and AF cou- 
plings for different form of initial randomness (Gaussian, 
symmetric and asymmetric rectangular) and for compar- 
ison calculations on the S = 1 model are also performed. 
The calculated distributions of the gaps are presented in 
Fig. EH 

As seen in Fig. ^| the slopes of the low-energy tail 
of the gap-distributions are approximately constant, and 
for our finite systems they are consistent with a vanishing 
gap exponent: 

uj^O (c! = 3). (25) 

During renormalization, as in 2c!, there is a large spin 
formation and the corresponding moment exponent is 
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FIG. 10: Probability distribution of the energy gap on the cu- 
bic lattice with mixed F and AF bonds, (a) Gaussian distri- 
bution, a = 1; (b) symmetric rectangular distribution (r = 0); 
(c) asymmetric rectangular distribution (r = 0.25); (d) S — 1 
symmetric rectangular distribution. The low-energy tails of 
the gap distributions for all cases, indicated by straight lines, 
have a slope -1, corresponding to uo = 0. 




FIG. 11: Scaling of the reduced gap distribution, P(L Z A) = 
i~ 2 Pt(A), for randomly frustrated 3d systems: a) Gaussian 
randomness, a = 1, b) uniform randomness. In both cases it 
is z = 1.5. 



C = 0.55, for symmetric distributions (Gaussian and rect- 
angular) and £ = 0.58 for the asymmetric rectangular 
distribution. Thus ( appears to be close to 1/2 in both 
cases. We have also studied the scaling behavior of the re- 
duced gap distribution, P(L Z A) = L~ 2 P L (A). In FiglTD 
we show a scaling collapse of the distributions, which is 
obtained by z ~ 1.5 independently of the disorder dis- 
tribution. The scaling curves seem to tend to a finite 
limiting value at A = 0, implying a gap exponent lu w 0. 
We can thus conclude that — within the range of valid- 
ity of the SDRG method — the relation in Eq. (|13[1 is not 
valid for frustrated 3d models. 



2. The J1-J2 model 

We also considered frustration caused by a competition 
between nearest and next-nearest neighbor couplings in 
order to see to what extent the universality of the SG 
phase, observed in 2d models, is valid in 3c?. Here we 
study systems at different points of the phase diagram, 
jmax^/jmax an( j f or different initial disorder D, using the 
same notations as for 2d. Typical gap distribution are 
shown in Fig. 1121 where we observe that the low-energy 
tail of the distributions in each case has approximately 
the same slope close to —1, which results in a gap ex- 
ponent, u) ss 0. This result is consistent with Ea. l|25() 
obtained for randomly frustrated models. During renor- 
malization large spin formation is observed, and the mo- 
ment exponent, £, is found to depend on the position in 
the phase diagram: for J™ ax /J™ ax = 0.5 and 1.0 it is 
given by C = -58 and .78, respectively (D = 1 is in both 
cases). The dynamical exponent in these cases was about 
the same (z » 3/2) as for randomly frustrated models. 
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FIG. 12: Probability distribution of the energy gap of the J\- 
J 2 model, (a) JJ ,ax /Jr ax = 0.5, D = 1; (b) J^/Jf^ = 1, 
D = 2; (c) J^ nax / J{ nax = 1, D = 1. (The slope of the straight 
lines in all cases is — 1.) 

Thus we can conclude that also in 3d the low-energy 
fixed points of random Heisenberg systems with differ- 
ent types of frustration are controlled by the same type 
of SGFP, having the same gap exponent, ui w 0, as in 
a two-dimensional SGFP. Therefore we conjecture that 
the ground states of these 3d frustrated models are in a 
spin glass phase, too. At these SGFP-s the dynamical 
exponent is constant, however the moment exponent has 
a system dependence. Thus the low-energy excitations 
have a universal scaling behavior, but the thermodynam- 
ical singularities in Ea. l|17|) . which depend on the value 
of C are system dependent. 



V. DISCUSSION 

In this paper we considered higher dimensional HAF-s 
and studied the effect of strong randomness on their low- 
energy/low-temperature properties by a numerical appli- 
cation of the SDRG method. Comparing with the known, 
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partially exact results for Id HAF-s we noticed several 
important differences. First, in higher dimensions one 
observes a strong universality scenario: there are only a 
few relevant fixed points (most important are the ran- 
dom AF fixed point and the SGFP) and their properties 
do not depend on the coordination number, the strength 
of disorder, value of the spin, etc., but just on the di- 
mension of the model and the degree of frustration in 
the system. In contrast to this, in random spin chains 
and ladders one has usually a continuum of low-energy 
fixed points parameterized by the value of the dynamical 
exponent z and which do depend on the aforementioned 
details. Second, in higher dimensional HAF-s the sin- 
gularities are controlled by (a few) conventional random 
fixed points, at which the dynamical exponent is finite. In 
higher dimensional systems there are no IRFP-s that can 
generally be found in (quasi) Id systems at random quan- 
tum critical points. A third difference between Id and 
higher-dimensional AF is the following: In Id the renor- 
malization of random spin-1/2 AF spin chains and ran- 
dom quantum ferromagnets, such as the random trans- 
verse Ising modelSi^i, lead to similar IRFPs at quantum 
critical points. In higher dimensional random transverse 
Ising models the random quantum critical point is still 
an IRFP^Si^, whereas for the random HAF, even at ran- 
dom quantum critical points, we found in this work the 
dynamical exponent to be finite. 

One remarkable aspect of our results is the observed 
universality of the fixed point controlling the low-energy 
characteristics of random frustrated systems 4 ^. The gap 
exponent of this so-called spin glass fixed point (SGFP) 
is numerically very close to zero— ^ and we can explain this 
observation in the following way. During rcnormalization 
there is a large spin formation in these systems and there- 
fore we expect that the low-energy excitations in d > 2 



are extended over the whole (finite) volume of the sys- 
tem (in a Id topology these excitations are not extended 
since unfavorable domains usually restrict the size of exci- 
tations). As a consequence these excitations can be con- 
sidered as compact objects so that their reduced (scale 
invariant) probability density P(L Z A) — L~ z Pl(A) has 
no size dependence for a fixed small gap, A. This last 
statement is consistent with a vanishing gap exponent, 
uj = 0, according to Eg. 111(1 and is supported by the 
numerical results in Fie llll 

Finally we make a remark about the accuracy of the 
results obtained by the SDRG method. It is generally 
expected that the SDRG method leads to asymptotically 
exact relations concerning singularities and scaling func- 
tions at IRFP-s 6 ' 41 . However, the same type of asymp- 
totic accuracy of the results is predicted at GFP-s and 
checked numerically by the density matrix renormaliza- 
tion group method^. Therefore we expect the predic- 
tions of the SDRG method about LSFP-s and the SGFP 
also to be correct. This expectation finds support in 
the results of numerical calculations for the Id Heisen- 
berg model with mixed F-AF couplings^ 6 - and for the ± J 
square lattice HAF—. Nevertheless alternative calcula- 
tions are necessary to check the validity of the predictions 
of our SDRG results. 
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